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Abstract 

Background: Scleractinian corals and their algal endosymbionts (genus Symbiodinium) exhibit distinct bathymetric 
distributions on coral reefs. Yet, few studies have assessed the evolutionary context of these ecological distributions 
by exploring the genetic diversity of closely related coral species and their associated Symbiodinium over large 
depth ranges. Here we assess the distribution and genetic diversity of five agariciid coral species {Agoricio humilis, 
A. agaricites, A. lamarcki, A. grahamae, and Heliosehs cucullata) and their algal endosymbionts {Symbiodinium) across 
a large depth gradient (2-60 m) covering shallow to mesophotic depths on a Caribbean reef. 

Results: The five agariciid species exhibited distinct depth distributions, and dominant Symbiodinium associations 
were found to be species-specific, with each of the agariciid species harbouring a distinct ITS2-DGGE profile (except 
for a shared profile between A. lamarcki and A. grahamae). Only A. lamarcki harboured different Symbiodinium types 
across its depth distribution (i.e. exhibited symbiont zonation). Phylogenetic analysis {atp6) of the coral hosts 
demonstrated a division of the Agaricia genus into two major lineages that correspond to their bathymetric 
distribution ("shallow": A. humilis I A. agaricites and "deep": A. lamarcki I A. grahamae), highlighting the role of 
depth-related factors in the diversification of these congeneric agariciid species. The divergence between "shallow" 
and "deep" host species was reflected in the relatedness of the associated Symbiodinium (with A. lamarcki and A. 
grahamae sharing an identical Symbiodinium profile, and A. humilis and A. agaricites harbouring a related ITS2 
sequence in their Symbiodinium profiles), corroborating the notion that brooding corals and their Symbiodinium are 
engaged in coevolutionary processes. 

Conclusions: Our findings support the hypothesis that the depth-related environmental gradient on reefs has 
played an important role in the diversification of the genus Agaricia and their associated Symbiodinium, resulting in 
a genetic segregation between coral host-symbiont communities at shallow and mesophotic depths. 
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Background 

Historically, coral species were thought to be extremely 
ecologically and morphologically plastic, because of the 
broad variation observed "within species" [1]. Advances 
in coral taxonomy and ecology have demonstrated that 
this variation often reflects sibling species occupying dis- 
tinct ecological niches [2]. The recent advent of molecu- 
lar approaches has exposed a further wealth of cryptic 
diversity in scleractinian corals [3-5] often in relation 
to particular environments, highlighting the role of 
niche partitioning [6-9] in addition to ecophenotypic 
plasticity [10,11]. Likewise, the algal endosymbionts 
that associate with reef-building corals and allow them 
to exploit light as an energy source in oligotrophic 
conditions, were originally thought to consist of a sin- 
gle algal symbiont, Symbiodinium microadriaticum 
[12]. By now, we know that Symbiodinium diversity in 
tropical scleractinian corals consists of several different 
phylogenetic clades [13-15], species [16,17] and hundreds 
of different subclades based on the ITS2 region of the 
ribosomal DNA [14,16]. These different Symbiodinium 
genotypes often represent physiologically distinct entities 
that can be strongly partitioned across host species, 
geographic regions and particular reef environments 
[18-23]. Although ecological niche partitioning across 
environmental gradients (e.g. light, temperature, turbi- 
dity) appears fundamental to maintaining the diversity 
of both scleractinian corals and their endosymbiotic di- 
noflagellates, the evolutionary context of this host- 
symbiont niche partitioning remains poorly understood. 

Despite the great diversity of Symbiodinium types 
found in reef invertebrates, many coral host species 
appear to associate with a specific genetic lineage (or 
group) of Symbiodinium [16]. Particularly in brooding 
corals, where the onset of symbiosis occurs through the 
vertical transmission of Symbiodinium from the mater- 
nal colony to the larvae, high levels of host-symbiont 
specificity have been observed [8,24-27]. Molecular stu- 
dies targeting both symbiotic partners have identified 
tightly coupled host-symbiont associations in brooding 
corals at the family [20], genus [9] and species level 
[8,9,26]. As such, natural selection may be acting on the 
level of the coral holobiont (i.e., host and symbiont com- 
bined) in species with a vertical symbiont acquisition 
mode, resulting in coevolution of both symbiotic part- 
ners. Incongruence between the phylogenies of corals 
and their associated Symbiodinium indicates that recom- 
bination of host-symbiont associations does occur over 
the course of evolution [13,16,28,29], however the occur- 
rences of such events may be less likely over ecological 
timescales. It is therefore important to study the evo- 
lutionary origin and implications of such host-symbiont 
specificity, particularly in light of the environmental 
change predicted over the next decades [30]. 



The evolutionary context of coral host-symbiont niche 
partitioning is also important in the context of reef 
connectivity between distinct environments and more 
specifically the ability of reef habitats that offer some 
kind of protection against certain disturbances, such as 
deep reefs, to act as a "protected" local source of 
propagules [31,32]. Although the association of a coral 
species with different Symbiodinium types may facili- 
tate its distribution across distinct reef environments 
[20,33-35], larval connectivity between these alternative 
environments may be hampered if symbiotic associa- 
tions are constrained through vertical symbiont trans- 
mission and recruits are not able to adopt symbiont 
types better-suited to their "new" environment [26]. 
Along those same lines, host populations may be lo- 
cally adapted to opposing environments and distinct 
symbiont populations may in fact reflect genetic diffe- 
rentiation of both symbiotic partners [31]. Such ge- 
netic partitioning of host and symbiont populations has 
been observed across shallow and deep populations of 
the pocilloporid coral species Madracis pharensis in the 
Caribbean [9] and Seriatopora hystrix on the Great 
Barrier Reef [8,26,32]. To gain a better understanding 
of the ability of mesophotic reefs to act as refugia (for 
species that occur over large depth ranges) and a subse- 
quent source of reproduction for degraded shallow reefs 
it is important to determine the overlap in genetic diver- 
sity of host and symbiont communities in shallow and 
deep reef environments. Given the prevalence of a 
brooding reproductive mode (>65%) and therefore verti- 
cal symbiont transmission on Caribbean reefs [31,36], 
there is a potential for strong genetic differentiation 
between host-symbiont associations in shallow and 
mesophotic communities in this region [27] . However, this 
hypothesis remains largely untested, primarily as molecu- 
lar data from mesophotic coral communities are limited. 

Mesophotic reef communities in the Caribbean are 
commonly dominated by members of the Agariciidae 
[37-39], and in particular by the species Agaricia 
lamarcki and Agaricia grahamae [31,40]. The plating 
growth forms, such as those of Agariciidae, are thought 
to be particularly well suited to maximize light capture 
in these low-light environments [41], Nonetheless, spe- 
cies of the genus Agaricia can be abundant throughout 
the entire depth range of Caribbean reefs and A. humilis, 
for example, occurs only in very shallow sections of reefs 
[42]. Species of Agaricia can exhibit high linear growth 
rates [43] and these brooding species are major con- 
tributors to the pool of coral recruits on Caribbean reefs 
[38,42,44,45]. Data on the phylogenetic affiliation of 
Symbiodinium associated with Western-Atlantic agari- 
ciids is limited, however records so far demonstrate 
association with Symbiodinium from clade C and occa- 
sionally clade D [16,46-50]. It remains unclear whether 
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the brooding reproductive mode of Western- Atlantic 
agariciids has resulted in host-symbiont specificity and 
to what extent the associated Symbiodinium diversity 
allows this genus to be abundant across large depth 
gradients. 

Here, we assessed Symbiodinium genetic diversity for 
the five most common agariciid corals in Curacao over a 
large depth range (2-60 m) encompassing both shallow 
and mesophotic depths. We addressed the level of 
specificity between the different coral species and their 
dominant Symbiodinium (i.e. the algal type(s) that are 
most abundantly present in the coral tissue), and eva- 
luated whether these host-symbiont associations are 
partitioned over depth. Additionally, sequencing of a 
host mitochondrial region disclosed agariciid evo- 
lutionary patterns, and allowed evaluation of further 
genetic structuring of the coral host in relation to sym- 
biont and depth. Results are discussed in light of exten- 
sive previous work on the brooding species of the coral 
genus Madracis [9,21,34,51] to address common 
patterns of genetic diversity over depth, and explore 
potential processes of codiversification in brooding 
Caribbean corals and their algal symbionts. Overall, the 
potential role of depth-related processes in the evolu- 
tion of brooding corals and their associated Symbio- 
dinium is discussed, as are the implications with 



regards to larval connectivity between shallow and deep 
coral reef habitats. 

Methods 

Study site and environmental conditions 

Collections were performed between 2008 and 2009 at 
the well-studied site Buoy Zero/One (Bak 1977; Bak 
et al. 2005; Vermeij et al. 2007; Frade et al. 2008a) in 
Curacao, southern Caribbean (12°07'3r N, 68°58'27"W) 
(Figure la). Information obtained from additional sam- 
ples of previous collections at the same study site in 
2004 and 2007 were added to the dataset. The study site 
is, as are many leeward sites on Curacao, characterized 
by a shallow terrace (50-100 m wide) that slopes down 
to 8-12 m, followed by a steep drop-off with a slope 
extending down to a deep terrace at 50-60 m depth [52]. 
Environmental conditions, including seawater tempera- 
ture and light attenuation were characterized during the 
beginning of the sampling period and described in Frade 
et al. [9]. Downwelling light availability (or PAR irra- 
diance) decreases exponentially with a K d s of 0.07 m _1 
[9], leading to -87% of surface irradiance at 2 m, 70% at 
5 m, 50% at 10 m, 17% at 25 m, 6% at 40 m, 3% at 50 m, 
and 1.5% at 60 m. Temperature at the study site shows 
seasonal variation, with mean seawater temperatures 
decreasing over depth (0.02°C per metre) [9]. The deeper 
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reef slope (40-60 m) experiences cold-water influxes due 
to thermocline movements, which can cause a sudden 
decrease in temperature of several degrees [9,39]. 

Agariciid species 

The five most common agariciid species in Curacao, 
Agaricia humilis Verrill, 1901, A. agaricites (Linnaeus, 
1758), A, lamarcki Milne Edwards and Haime, 1851, A. 
grahamae Wells, 1973, and Helioseris cucullata (Ellis 
and Solander, 1786) were included in this study. Of the 
other three Caribbean agariciid species, A undata (Ellis 
and Solander, 1786) is rare, and A, fragilis Dana, 1848 
and A, tenuifola Dana, 1848 are absent on the reefs of 
Curacao. Coral species were identified following the 
taxonomic features specified by Wells [53], Veron and 
Stafford-Smith [54] and Humann & Deloach [55] 
(Table 1). Fragments of colonies that were collected but 
that did not clearly conform to these species descrip- 
tions were labeled as "Agaricia sp.". 

Species abundances and sampling approach 

Species abundances of the different agariciid species 
were determined in 2009 by counting all individuals in 
30-50 quadrats (1 x 1 m) laid out randomly (i.e. using a 
list of randomly generated numbers) along a 100 m long 
transect parallel to the shore at each sample-depth (2, 5, 
10, 25, 40, 50, and 60 m). A total of 335 agariciid coral 
colonies were sampled over a depth range from 2 to 
60 m from a single reef location (Buoy Zero/One), by 
removing a small fragment from the colony (~3 cm 2 ). 
Fragments were collected at depth intervals that reflect 
the bathymetric distributions of each species: A, humilis 
at 2, 5 and 10 m (n = 71); A. agaricites at 5, 10, 15, 25, 30, 
40 and 50 m (n = 113); A. lamarcki at 10, 25, 40 and 
50 m (n = 105); A. grahamae at 50 and 60 m (n = 20); and 
H. cucullata at 25 and 40 m (n = 16). Tissue from col- 
lected coral samples was separated from the skeleton 



using a modified airgun attached to a SCUBA cylinder or 
with a sterilized scalpel, and was subsequently stored in 
20% DMSO/EDTA preservation buffer or 95% ethanol at 
-20°C until further processing. DNA was extracted 
from the tissue using a Qiagen Plant Mini Kit, MoBio 
Ultra Clean Soil DNA Kit (following the manufacturer s 
instructions), or a slightly modified method used for 
black tiger shrimp [32,56]. 

Symbiodinium ITS2 analyses 

The internal transcribed spacer (ITS2) region of 
Symbiodinium rDNA was amplified for all specimens 
(n = 335) using Symbiodinium-specific primers [46] as 
described in Bongaerts et al. [8]. The amplified ITS2 
fragments were separated using denaturing gradient gel 
electrophoresis (DGGE) on a CBScientific System 
following conditions outlined in Sampayo et al. [20], to 
identify the dominant Symbiodinium types in each 
sample (ITS2-DGGE fails to detect symbionts present 
in low abundances [57,58]). On each DGGE gel, we ran 
representative samples to allow for profile cross- 
comparison and identification. Representative, domin- 
ant bands of each characteristic profile were excised 
(from several replicate profiles), eluted overnight in 
dH 2 0, re-amplified and purified (using ExoSAP-IT) 
prior to sequencing. The re-amplified PCR products 
were sequenced in both the forward and reverse direc- 
tions (ABI BigDye Terminator chemistry, Australian 
Genome Research Facility). All chromatograms were 
analyzed using Codoncode Aligner, with sequences 
being aligned with MUSCLE and blasted on GenBank 
(http://www.ncbi.nlm.nih.gov/BLAST). Phylogenetic ana- 
lyses of sequences were performed using maximum parsi- 
mony and maximum likelihood in MEGA 4 [59] under 
the delayed transition setting and calculation of bootstrap 
support values based on 1000 replicates. 



Table 1 Morphological characteristics of the five studied Agariciid species 



Species 


n 


Colony 
shape 


Plates 


Cora Mites 


Valley orientation 


Septocostae 


Other 


Agaricia humilis 


71 


Submassive, 
encrusing, 
or plating 


n.a. 


Deep, monocentric 

corallites (or 2-3 centres together) 


Centres separated by collines 
of irregular height running in 
irregular directions 


Alternating 




Agaricia agaricites 


113 


Encrusting, 
or plating 


Bifacial 
(fronds) 


Short, deep valleys 
(usually <50 mm). 


Irregular, or concentric on 
plates 






Agaricia lamarcki 


105 


Plating 


Unifacial 


Valleys with widely 
spaced centres 


Concentric valleys 


Alternating 


White & star- 
shaped mouths 


Agaricia grahamae 


20 


Plating 


Unifacial 


Long, V-shaped valleys 


Roughly concentric valleys 


Non- 
alternating 




Helioseris cucullata 


16 


Plating 


Unifacial 


Non-continuous valleys 
with centres nested against 
steep ridges 


Roughly concentric valleys 


Alternating 





Taxonomic features of the different Caribbean Agariciid species studied here, as specified by Wells [53], Veron and Stafford-Smith [54] and Humann & 
Deloach [55]. 
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Host atp6 sequence analyses 

The mitochondrial atp6 gene (and parts of the flanking 
nad6 and nad4) was amplified for a subset of samples 
(n = 93: ~5 samples were randomly picked per species 
per depth) using the newly designed primers Aga-atp6- 
Fl: GGCTTTATTTGGGGCTGAA and Aga-atp6-Rl: 
CCC AC AAAACC AAAG C ACTTA. Primers were desig- 
ned from the mitogenomes of A. humilis (NC008160) 
and Pavona clavus (NC008160) [60]. PCR amplifications 
were performed with 1.0 ul of DNA, 2.0 ul lOx PCR 
buffer (Invitrogen), 1.0 ul 50 mM MgCl 2 , 0.7 ul 10 mM 
dNTPs, 0.7 ul 10 mM Aga-atp6-Fl, 0.7 ul 10 mM Aga- 
atp6-Rl, 0.15 ul of Platinum Taq DNA Polymerase 
(Invitrogen) and dH 2 0 water to a total volume of 20 ml 
per reaction. The cycling protocol was: 1 x 95°C (5 min); 
30 x [95°C (30 s), 58°C (60 s), 70°C (90 s)]; 1 x 70°C 
(10 min). PCR products were purified, sequenced, 
aligned and analysed as specified for the Symbiodinium 
ITS2 sequences. Additional sequences were retrieved 
from GenBank to compare the atp6-based phylogenies 
with those based on the mitochondrial cytochrome 
oxidase 1 (COI) and the nuclear 28S region [60-63]. 

Phylogenetic analyses of sequences were constructed 
using maximum likelihood (ML), maximum parsimony 
(MP) and Bayesian (BAY) methods through the pro- 
grams MEGA5 (Tamura et al. 2011) and MrBayes 
(Huelsenbeck and Ronquist 2001) respectively. The 
best-fit model of molecular evolution was selected by 
hierarchical Akaike information criterion (AIC) using 
jModeltest [64] with a GTR model with invariant sites 
best describing the atp6 data under a log likelihood 
optimality criterion. For the COI and 28S sequences 
retrieved from GenBank, a GTR model with gamma- 
shaped variation and invariant sites best described the 
data of both these regions. Maximum likelihood ana- 
lyses were performed using 1000 non-parametric boot- 
strap replicates. Bayesian analyses were performed with 
the Markov Chain Monte Carlo search run with 4 
chains for 10 6 generations, with a sample frequency of 
100 generations and a "burn-in" of 2500 trees. 

Statistics 

The contributions of taxonomic species, depth and 
Symbiodinium type on host genotypic variability (atp6) 
were assessed for positively identified Agaricia specimens 
under the AMOVA framework using GenALEx V5 [65], 
with either depth or Symbiodinium type nested within 
taxonomic species. In host species for which multiple 
Symbiodinium profiles were observed, an (nested) analysis 
of similarity (ANOSIM) was carried out using Bray-Curtis 
Distance in the software package Primer v6 to test for 
differences between sampling years (Two-way ANOSIM 
with depth nested within year) and depths (One- Way 
ANOSIM comparing individual depth groups). 



Results 

Species abundances over depth 

The species abundances at Buoy Zero/One differed 
among depths and revealed different distribution ranges 
for the five agariciid species (Figure lc). At the shallowest 
depths (2-5 m), A. humilis is the most dominant species, 
with A. agaricites starting to occur at 5 m. At 10-15 m, A. 
agaricites becomes the most dominant species, with 
colonies of A. humilis and A. lamarcki being present in 
lower abundances. At 25 m depth A. agaricites and A. 
lamarcki occur in roughly equal abundances, but A. 
lamarcki takes over as the dominant species beyond this 
depth. From 40 m onwards, A. grahamae starts occurring, 
although in low numbers. H. cucullata was observed from 
10-60 m, but always in relatively lower abundances. 

Symbiont diversity across host species and depth 

Across the five different host species, a total of six dis- 
tinct ITS2-DGGE profiles were distinguished (Figure 2; 
see Additional file 1 for DGGE profiles). A distinct pro- 
file was observed for A. humilis (PI, n = 71), A. agari- 
cites (P2, n = 113), and H. cucullata (P5, n = 16). Three 
different profiles (P3, n = 34; P4, n = 69; P4*, n = 2) were 
observed for A. lamarcki, of which one (P4) was shared 
with the species A. grahamae (n = 20) (Figure 2). Given 
that five of the six profiles contained at least three co- 
dominant ITS2 sequences, we decided to assign a num- 
ber to the DGGE profiles (PI, P2, P3, P4, P4* and P5) 
rather than referring to individual Symbiodinium types. 
Although these co-dominant sequences likely represent 
intra-genomic variants within the rDNA of a single sym- 
biont lineage (little variation in relative band intensity 
was observed) [66], it cannot be excluded that some of 
the profiles may represent a mix of distinct Symbio- 
dinium types. Additionally, other Symbiodium types may 
be present at background levels that cannot be detected 
due to the limitations of ITS2-DGGE [57,58]. 

All the ITS2-DGGE fingerprint profiles for Agaricia 
spp. contained the Symbiodinium C3 sequence, in com- 
bination with 2-3 sequences that were all closely related 
to C3 (Figure 3). All of these sequences were novel 
(GenBank Accession Numbers KF551185-KF551192), 
except for C3b in A. humilis, C3d in A. lamarcki, and 
Cll in A. lamarcki I grahamae. Novel sequences are 
identified (in this manuscript) by a nomenclature that 
specifies the sequence to which they are most related, 
followed by a capital N (indicating novel sequence) and 
an arbitrary number in italics (e.g. C3bA/7). Except for 
the shared C3 sequence, all profiles contained distinct 
Symbiodinium sequences (i.e. no overlap between 
profiles) indicating that each of these profiles contain (a) 
different Symbiodinium type(s). However, one Symbio- 
dinium ITS2 profile (P4*) was observed in only 2 
colonies of A. lamarcki (at 25 m) and was identical to 
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A. humilis 



A. agaricites 



A. lamarcki 



A. grahamae 



H. cucullata 
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5 m 
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n = 23 


10 m 

© 

n = 25 
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• 
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• 
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• 
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© 
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© 

n = 25 
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© 

n = 20 


© 

n = 15 


• 


• 


• 
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n = 32 


• 

n = 22 


• 

n = 11 


• 


• 


• 


• 


• 




• 

n = 11 


• 

n = 9 


• 


• 


• 


© 

n = 10 


© 

n = 6 


• 


• 



□ P1 

□ P2 

□ P3 
■ P4 

□ P5 



Figure 2 Distribution of Symbiodinium ITS2 profiles across the five agariciid host species and depths. Each pie chart represents the 
sampled population of a host species at a certain depth. Asterisk (*) indicates a P4 profile that had an additional CI band. Figure legend text. 



Symbiont - ITS2 

C3A/3 
A. agaricites 4>p 
C3A/2 



/~r>iC3N7 
QC3N6 




A. lamarcki 



C11A/4 

A. grahamae 
A. lamarcki 



QCN8 
H. cucullata 



Figure 3 Sequence network of Symbiodinium ITS2 types 
observed for each DGGE profile obtained from agariciid 
species. Ancestral Symbiodinium types CI and C3 are represented 
by squares. Circle size is not representative of frequency. Dashed 
lines indicate sequence gaps (with gap size in number of base pairs). 
Colours group the different ITS2 sequences observed in each DGGE 
profile. The C3 type (indicated by black square) is found in all DGGE 
profiles except for P5. 



the P4 profile but with an extra band that corresponded 
with the CI sequence. The several Agaricia morphs that 
could not be unambiguously identified {Agaricia sp.) 
all contained the P4 profile (n = 10). For H. cucullata 
we could only obtain a single recoverable ITS2 se- 
quence (CN8), which represented a more distant, 
novel Symbiodinium type (GenBank Accession Num- 
ber KF551192) (Figure 3). 

A humilis, A. agaricites, A. grahamae and H. cucul- 
lata all harboured a single symbiont profile across the 
entire depth range across which they were sampled. A. 
lamarcki, the only species observed to harbour multiple 
symbiont profiles, showed a marked zonation, with 
colonies at 10 and 25 m depth harbouring either P3 
(n = 34) or P4 (n = 36), whereas colonies (>40 m) exclu- 
sively harboured P4 (n = 33). Colonies hosting a mix of 
P3 and P4 profiles were not observed. There were no 
significant differences (Two-Way ANOSIM; depth 
nested within sampling year) between sampling years. 
The pairwise comparisons between the different depths 
confirmed that the symbiont community associated with 
A. lamarcki at 40 m was significantly different from 
10 and 25 m (One-Way ANOSIM; not taking into ac- 
count sampling year; respectively R = 0.15; p < 0.01 
and R = 0.082; p < 0.05). 

Host genetic structure across taxonomic species and 
depth 

Seven different mtDNA haplotypes were identified 
across the 93 specimens for which the atp6 region 
(1260 bp including indels) was sequenced (Figure 4; 
GenBank Accession Numbers KF551193-KF551199). A 
single haplotype was observed for H. cucullata (n = 10), 
and indels were only observed when Agaricia sequences 
were aligned with those of H. cucullata. The remaining 
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Figure 4 Unrooted network of host atp6 mitochondrial 
haplotypes. Circle size indicates the relative frequency of each 
haplotype. Pie chart indicates for which agariciid species the 
haplotype was observed (Ahum = A humilis (n = 16), Aaga = A 
ogaricites (n = 22), Alam = A lamarcki (n = 27), Agra = A grahamaei 
(n = 9), Asp. = Unidentified Agoricio species (n = 9), Hcuc = H. 
cuculloto (n = 10)), and colour indicates the Symbiodinium profile 
that was hosted. 



haplotypes were observed for Agaricia specimens and 
the number of substitutions between these sequences 
ranged from 0 to 16. The Agaricia haplotype network 
consists of two clades and four main haplotypes (i.e. 
subclades), with the majority of A, lamarcki and all A. 
grahamae specimens belonging to clade 2 (respectively 
90 and 100%), most A. agaricites belonging to subclade 
la (86%), and A. humilis belonging to subclades lb and 
lc (69%). The Agaricia sp. specimens that were se- 
quenced (n = 9 of 10) were found across the two main 
clades. Under the AMOVA framework (including only 
identified Agaricia specimens), 60% of molecular vari- 
ance (® S pp_ TO T = 0-599; P < 0.001) was explained by 
morphotaxonomic species designation. No significant 
contributions of Symbiodinium type (2%; ®sym-spp = 
0.058; P = 0.232) and depth (0%; O dep . S pp = -0.059; P = 
0.668) could be detected, when nested within taxonomic 
species. 

Phylogenetic analyses (ML, MP and BAY) of the atp6 
region of the Agaricia specimens (Figure 5), using H. 
cucullata and Pavona clavus as outgroups, supported 
the monophyly of the genus Agaricia and subdivision of 
the genus into two main lineages (with bootstrap values 



across phylogenetic methods higher than 97), with 
sequences of the two lineages largely corresponding with 
A. lamarcki I A. grahamae and A. humilis I A. agaricites 
specimens respectively. The latter lineage is further 
subdivided into two lineages largely corresponding with 
A. humilis and A. agaricites (ML/MP/BAY: 94/69/93). 
Despite the mostly invariant host sequences among con- 
specifics, a small number of specimens belonging to dif- 
ferent species were observed within each lineage. 
However, in all these cases the symbiont type hosted 
by these specimens confirmed the taxonomic identity 
(see Figure 4). Specimens that could not be identified 
{Agaricia sp.) were found to belong to the three diffe- 
rent groups, although all hosting Symbiodinium profile 
P4, with 5 specimens from 25-40 m in the A. lamarckil 
grahamae group, 2 specimens from 50 m in the A 
agaricites group, and 2 specimens from 60 m in the A. 
humilis group. 

Analyses of COI and 28S sequences for the same 
Agaricia spp. [60-63], using Pavona spp. and respectively 
Leptoseris sp. and H. cucullata as outgroups, corroborate 
the phylogenetic pattern revealed through the atp6 
region (Figure 5b, c). However, where COI shows no 
differentiation between A. humilis and A, agaricites, the 
28S region differentiates between A. lamarcki and A. 
grahamae specimens. As such, the 28S region supports 
the subdivision of the Agaricia genus into two main line- 
ages, corresponding to shallow (A. humilis and A. agari- 
cites) and deeper species (A. lamarcki and A. grahamae). 

Despite this cladal division into shallow and deep 
Agaricia species, no clear genetic subdivision was 
observed within species over depth (Figure 5), or in the 
case of A. lamarcki in relation to symbiont type 
(Figure 4). Nonetheless, several shallow specimens of A. 
lamarcki were found within the A. humilis (n = 3; 10 m) 
and A. agaricites groups (n = 1; 25 m). Additionally, 
several specimens of A. agaricites from 25 (n = 1) and 
50 m (n = 2) depth were observed in the A. lamarckil 
grahamae lineage. 

Discussion 

This study assessed the depth distribution and genetic 
variation of five agariciid coral species and their asso- 
ciated Symbiodinium across a large depth gradient 
(2-60 m) on a single reef. Overall, we found a tight 
coupling between coral species and their dominant 
Symbiodinium at the study location, and strong evi- 
dence for depth-related niche partitioning of these 
associations. Depth-related parameters appear to have 
played a major role in the diversification of the coral 
genus Agaricia and its dominant algal endosymbionts, 
resulting in a clear genetic partitioning between coral 
host-symbiont communities across shallow and meso- 
photic depths. 
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Figure 5 Phylogenetic trees of Agaricia spp. based on mitochondrial and nuclear markers. Phylogenetic trees (maximum parsimony) of 
Agoricio spp. based on (a) mitochondrial atp6 (this study), and (b) nuclear 28S and (c) mitochondrial COI (previous studies: Medina et al. 2006; 
Fukami et al. 2008; Shearer and Coffroth 2008; Barbeitos et al 2010). Povono spp., Helioseris cucullata and Leptoseris sp. are used as outgroups. 
Bootstrap values are based on Bayesian (BAY), maximum parsimony (MP) and maximum likelihood (ML) respectively, with only probabilities over 
50% shown. Depth range and number of specimens are mentioned for each species, coloured boxes indicates the Symbiodinium profile that was 
hosted. Sequences retrieved from GenBank are indicated by their Accession Number. 



Symbiodinium diversity and host-symbiont specificity 

Associations between the coral host and its dominant 
algal endosymbionts {Symbiodinium) at the study loca- 
tion were found to be specific in that each of the five 
agariciid species harboured a distinct Symbiodinium 
profile (except for A lamarcki and A. grahamae that 
had one profile in common) (Figure 2). The four Sym- 
biodinium profiles found in association with the four 
Agaricia species all shared the C3 sequence as an 
intragenomic variant and contained closely related 
Symbiodinium sequences (Figure 3). This corroborates 
observations in other parts of the Caribbean (Bahamas, 
Barbados, Belize, Florida, Mexico, US Virgin Islands), 
where Agaricia species have been found in association 
with C3-related Symbiodinium types [16,46-50]. 

Nonetheless, except for the C3b-related Symbiodi- 
nium type in A. humilis (profile PI), the agariciids stu- 
died here hosted distinct Symbiodinium types compared 
to those previously reported (e.g. Bl, C3, C3a, C3q, 
Dla) [16,46-50]. A, lamarcki associated with two pro- 
files (P3 and P4) containing respectively a C3d and Cll 
sequence, which indicate Symbiodinium types related to 
those found in respectively Montastraea (C3d) and 
mussid and faviid genera such as Scolymia, Mussa and 
Mycetophyllia (Cll) [16,35,46]. H. cucullata hosted a 
Symbiodinium type (CA/S/profile P5) that is very distinct 



from the C3 type previously reported in association with 
this species in other parts of the Caribbean [46,48,50], and 
does not contain the C3 sequence as an intragenomic 
variant. Some of these differences in Symbiodinium associ- 
ated with Agaricia spp. (other studies) can be contributed 
to geographic separation of populations from the Western 
and Eastern Caribbean [50], however it also demonstrates 
that differences in symbiont associations exist even at 
smaller geographic scales, such as within the Eastern 
Caribbean (i.e., between Curacao and Barbados). These 
observations corroborate the notion that the level of host- 
symbiont specificity can vary across geographic locations 
[22,29,32], and therefore has to be assessed locally in 
terms of the ecological implications. 

The maternal transmission of endosymbionts leads to 
relative isolation of host-specific Symbiodinium popu- 
lations, and is likely to be responsible for the local host 
specificity and the geographic variation in host-symbiont 
associations, as Symbiodinium types coevolve locally 
with their coral partners [22]. Similar host-symbiont 
specificity has been observed in other coral genera with a 
vertical symbiont acquisition mode, such as in various 
pocilliporid genera (e.g. Seriatopora, Pocillopora, Stylo- 
phora, Madracis) and the genera Montipora and Porites 
[20,22-25,29]. Although the incongruence between host 
and symbiont phylogenies points towards the importance 
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of occasional host-symbiont recombination [13,16,28,29], 
such recombination through "symbiont replacement" is 
unlikely to occur over ecological timescales [67]. There- 
fore, despite the variation in host-symbiont associations 
observed for Agariciids across the broader Caribbean, it 
remains questionable whether the agariciid species at the 
study location can respond to rapid environmental change 
due to the lack of standing variation (i.e. diversity) in 
Symbiodinium associations. 

The host-symbiont specificity observed here relates to 
the dominant Symbiodium types that are hosted by each 
species, and high-resolution approaches to the detection 
of Symbiodinium (e.g. using quantitative PCR or next- 
generation sequencing) may demonstrate that these 
agariciid species also host non-specific background 
Symbiodinium [57,58]. In addition, we observed the ge- 
neralist Symbiodinium type CI in 2 out of 105 sampled 
specimens of A lamarcki (in addition to profile P4), 
indicating that non-specific associations occasionally do 
occur [29]. Although the potential of "symbiont shuf- 
fling" (i.e. shift in dominant Symbiodinium within a coral 
colony over its lifespan) as a response to environmental 
change remains debated [67], unusual associations (such 
as the CI hosted in A. lamarcki) may become more 
abundant if these exhibit higher levels of holobiont 
fitness under the predicted changing environmental 
conditions. 

Depth-partitioning of agariciid host species and 
associated Symbiodinium 

The five agariciid coral species exhibited a clear 
partitioning across the depth-mediated environmental 
gradient, progressing from respectively A humilis, A. 
agaricites and A. lamarcki as most abundant species 
from shallow to deep (Figure lc). A. humilis is a clear 
"shallow-water specialist" largely restricted to the shal- 
low reef terrace (2-10 m), whereas A. grahamae is a 
"deep-water specialist" that exclusively occurs in the 
mesophotic zone (40-60 m). In contrast, both A. agari- 
cites and A. lamarcki are more generalistic in terms of 
their depth distribution, albeit that A, agaricites predom- 
inates at shallower depths (5-25 m) and A. lamarcki at 
greater depths (25-60 m), with 25 m representing a tran- 
sitory depth between both species (Figure lc). These 
depth-distributions are from a single study location, and 
may vary across locations with different environmental 
conditions (e.g. water clarity) [68]. However, similar 
zonation of Agaricia species has been reported across 
the Caribbean [37,38,69-72], with A. humilis and A. 
agaricites most abundant on the shallow end of the 
depth spectrum, and A. lamarcki and A. grahamae 
predominating at greater depths. H. cucullata was found 
across a large depth range (10-60 m), however only in 
very low abundances, reflecting the remarkable decline 



of this species on the reefs in Curacao compared 
to -30 years ago [44,45]. 

Four of the agariciid species hosted a single 
Symbiodinium type over their entire depth range. In 
contrast, the depth-generalist species A. lamarcki 
exhibited symbiont depth zonation, changing from 
hosting P3 (i.e. C3/C3d -related) or P4 (i.e. C3/Cll-re- 
lated) in the shallow to exclusively P4 at mesophotic 
depths. The "deep-water specialist" A. grahamae har- 
boured this same P4 symbiont profile, indicating that 
this Symbiodinium type may be particularly well-adapted 
to deeper water conditions. Similar putative "deep-spe- 
cialist" Symbiodinium types have been identified in the 
genera Madracis and Montastraea [21,34,35]. Con- 
versely, the ability of A. lamarcki to associate with 
another Symbiodinium lineage (represented by the P3 
profile) may facilitate its broader depth distribution (i.e., 
extending into the shallow) compared to A. grahamae 
(although the P4 profile was found at shallow depths in 
A. lamarcki as well). Symbiont zonation appears to be a 
common feature of depth-generalist coral species in the 
Caribbean [31], however the occurrence of a single 
Symbiodinium profile in A. agaricites from 5-50 m depth 
indicates that such zonation, however, is not a prerequi- 
site for broad distribution patterns [27,73]. 

When considering both the average densities of host 
species and dominant symbiont associations across 
depths, it is obvious that host-Symbiodinium assem- 
blages exhibit a clear niche partitioning over depth, for 
all assemblages apart from H. cucullata I P5 (conceptual 
diagram shown in Figure 6). Although both A. humilis I 
PI and A. agaricites I V2 predominate in the shallow reef 
environment, A. humilis I VI thrives in the shallowest of 
depths (2-10 m) whereas A, agaricites I V2 is most abun- 
dant at intermediate depths (10-25 m). As such, A 
humilis I VI is restricted to depths with high-irradiance 
conditions (87-50% of surface irradiance) and A. aga- 
ricites I V2 occurs in all but the brightest light conditions 
(1.5-70%). At greater depths (> 25 m), however, A. 
lamarcki/V^ becomes the most dominant host-symbiont 
assemblage. Given that all sampled colonies of A. 
lamarcki at 10 m depth (n = 40) were growing in cryptic 
locations on the reef (Bongaerts and Frade pers. obs.), 
the shallow end of A. lamarcki s depth distribution is 
clearly determined by a maximum tolerable irradiance. 
Additionally, the A. lamarcki/V3 host-symbiont asso- 
ciation, which only occurs at 25 m and cryptically at 
10 m, occupies a niche that encompasses a very narrow 
light range. Instead, A. grahamae that shares symbiont 
profile P4 (C3/Cll-related) with A. lamarcki, starts to 
appear at slightly greater depths than the A. lamarcki/V4 
combination, indicating the potential (acclimatization) 
importance of the host component in mediating the 
distribution of this host-symbiont partnership [34]. 
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Although A. lamarcki was found to be most abundant at 
the deeper end of the sampled depth spectrum (60 m), 
anecdotal observations at greater depths (65-85 m) 
indicate that A. grahamae may in fact become the more 
dominant species (Figure 7; Vermeij and Bongaerts, 
pers. obs.). 

Although incident irradiance is likely an important fac- 
tor to explain the observed depth partitioning of host- 
symbiont associations, other environmental variables, 
such as light spectral composition [74], temperature [75] 
and nutrient/plankton availability [26] may be important 
as well For example, the depth range between 25 and 
40 m appears to be a transition zone of strong genetic 
change, which corresponds to the depth range where 
influxes of deep oceanic water become noticeable in the 
temperature regime [21,39]. Given the various abiotic 
factors that vary over depth it remains impossible to 
evaluate the individual contributions of these factors in 
the observed partitioning, however light and tempe- 
rature are likely to play important roles given their 
dramatic changes over depth at the study site [21,76]. 




Figure 7 Photo of Agaricia community beyond the depth range 

of this study. Paper-thin Agaricia grahamae communities at 86 m 

depth on the leeward shore of Curacao. 
\ ) 



While representing a massive undertaking, repeating 
similar symbiont genotyping efforts across a range of 
different locations would provide us with further insight 
into the role of different environmental conditions in the 
depth partitioning of host-symbiont associations. 

Host and symbiont (co) evolution 

The mtDNA-based phylogenetic analysis showed that 
none of the Agaricia species studied is monophyletic, how- 
ever, the presently accepted species designation can 
explain the majority of the observed molecular variance 
(-60%). Two major lineages were observed for the 
Agaricia specimens based on the atp6 region, basically di- 
viding the genus into a shallow {A. humilis/A. agaricites - 
Clade 1) and a deep (A. lamarckilA. grahamae - Clade 2) 
species group. The atp6 phylogeny is consistent with the 
phylogeny of the ribosomal 28S region (outgroup H. 
cucullata not included) produced with GenBank se- 
quences, although the 28S region appears more variable 
(Figure 5). Although the atp6 region does not separate A. 
lamarcki and A. grahamae, the phylogeny of the 28S 
region demonstrates that these are unlikely to represent 
synonyms. The H. cucullata clade is genetically very 
distinct and even shows (across markers) a closer relation 
to the Indo-Pacific genus Pavona than the Atlantic 
Agaricia genus (Figure 5). The phylogeny of Agaricia, 
forming two major clusters representing shallow and deep 
species, indicates that depth-related factors have played a 
major role in the original divergence within this genus. 
Given that also within each of these clusters, each species 
again exhibit distinct depth distributions on respectively 
the shallow and deep reef (Figure lc), highlights the poten- 
tial role of depth-related factors in the diversification 
within this genus throughout its evolutionary history. 

Within the genus Agaricia all of the four major host 
atp6 haplotypes are, in fact, shared by more than one 
Agaricia taxon. This shared genetic polymorphism and 
repetitive non-monophyly has been shown in the past 
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for many other coral taxa [9,77-81]. Such incongruence 
between morphospecies and genetically delimited phy- 
togenies have been explained by several authors as being 
the result of phenomena such as morphological conver- 
gence (homoplasy), phenotypic plasticity, recent speci- 
ation with incomplete lineage sorting, morphological 
stasis, and/or interspecific introgressive hybridization 
[80]. Although introgressive hybridisation has been sug- 
gested in other brooding coral species [9,80,82], the fact 
that we only used a single mitochondrial marker 
prevents inference of such processes in the current 
study. The "misplaced specimens" in the phylogeny all 
hosted the Symbiodinium profile specific to that species 
(Figure 2), confirming that these are unlikely to repre- 
sent taxonomic misidentifications. However, the obser- 
vation that several shallower specimens of A lamarcki 
(n = 4) grouped with the "shallow" species group (Clade 
la-c), and several deeper specimens of A. agaricites (n = 3) 
grouped with the "deep" species group (Clade 2) 
(Figure 5), could potentially be indicative of hybridization 
due to low densities of conspecific sperm. 

The evolutionary divergence between "shallow" and 
"deep" Agaricia host species and genetic relatedness 
among them is reflected in the symbiont associations, 
with A. lamarcki and A. grahamae sharing an identical 
Symbiodinium profile (P4), and A. humilis and A. agari- 
cites harbouring a related ITS2 sequence in their 
Symbiodinium profiles (C3b in PI and C3A/2 in P2) 
(Figure 3). This highlights, despite potential hybridisa- 
tion within the Agaricia genus (which could result in the 
crossing over of Symbiodinium types), the diversification 
of both host and symbiont through reciprocal evolution- 
ary changes (at least at a microevolutionary scale). These 
observations corroborate the growing body of evidence 
suggesting that corals with a vertical symbiont acquisi- 
tion (i.e., usually brooding corals) and their associated, 
dominant Symbiodinium are involved in revolutionary 
processes [8,9,26,29,83]. 

A common pattern for Caribbean brooding corals? 

The depth niche partitioning, host phylogenetic struc- 
ture and local host-symbiont specificity observed for the 
genus Agaricia has strong similarities to that of the 
model genus Madracis, which has been extensively stud- 
ied at the exact same study site [9,21,34,51,76,83-87]. 
The Caribbean genus Madracis comprises six common 
taxonomic species whose distributions are partitioned 
over the reef slope, including shallow (M. mirabilis and 
M. decactis) and deep (M. carmabi and M. formosa) spe- 
cialists and two depth-generalist species (M. pharensis 
and M. senaria). However, similar to A. lamarcki, the 
distribution of M. pharensis and M. senaria is restricted 
to cryptic locations at depths of 10 m and shallower 
[51]. As such, the genus Madracis can be roughly 



divided into "shallow" species (M. mirabilis and M. 
decactis) and "deep" (M. pharensis, M. senaria, M. 
carmabi and M. formosa) species. This subdivision cor- 
responds to the two main clades observed in the 
Madracis phylogeny based on the mitochondrial nadS 
region [9], similar to the genetic divergence observed 
between "shallow" and "deep" species in the genus 
Agaricia, 

With regards to the associated Symbiodinium, the 
genus Madracis harbours three physiologically distinct 
but phylogenetically closely related symbiont types (B7, 
B13, B15), which are also partitioned over depth and to 
a certain extent specific to certain host species [21,34]. 
The shift in symbionts observed for M. pharensis corre- 
sponds to distinct physiological capacities of the symbi- 
onts involved [34], but also corresponds with a genetic 
divergence between shallow and deep populations of the 
host species M. pharensis [9]. Such intra-specific genetic 
divergence was not observed for A. lamarcki (the only 
Agaricia species that exhibited symbiont zonation), how- 
ever it does corroborate the general pattern of strong 
host-symbiont coupling in Agaricia (this study). 

Overall, the observed niche partitioning and evolution- 
ary split between shallow and deep host lineages 
observed for the genera Agaricia and Madracis, high- 
light the role of the depth-related environmental gradi- 
ent (i.e. differences in light, temperature, etc.) as a 
driving factor in the evolution of the congeneric species 
associated with these genera. The molecular phylogenies 
of both genera are relatively incongruent relative to the 
taxonomic classification based on morphological traits 
[9,84], highlighting the potential additional importance 
of introgressive hybridization. However, the role of intro- 
gressive hybridization in Agaricia is hard to infer given 
the use of a single, maternally-inherited marker region 
in this study. Nonetheless, the fact that species in both 
genera remain recognizable and have distinct ecological 
distributions does point towards the importance of 
divergent/disruptive selection across a depth gradient in 
the diversification of these genera. Given the prevalence 
of a brooding reproductive strategy [36] and vertical 
symbiont acquisition mode in the Caribbean, ecological 
adaptation to different ends of the depth spectrum may 
have played a major role in the diversification of many 
Caribbean corals. 

Implications for vertical connectivity 

In recent years, there has been a growing interest in the 
role of deep reef areas as local refugia and a source of 
propagules [31], particularly given the observed relative 
stablity of some mesophotic reef environments [39] and 
accumulating evidence of localized recruitment on coral 
reefs [88-90]. The extent of genetic specialization within 
the genus Agaricia to both shallow and deep reef 
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environments (Figure 6) indicates that, although the 
mesophotic Agaricia communities have remained rela- 
tively stable in the geographic location of our study [39], 
these communities are unlikely to provide significant 
amounts of propagules to ensure rapid shallow-reef 
recovery post-disturbance. The unrelated genus Madracis, 
shows a similar genetic segregation between its host- 
symbiont communities at mesophotic (dominated by M. 
pharensis I B15 and M.formosa I B15) and shallow depths 
[9,21], confirming that this pattern may extent across a 
broad range of coral genera. Nonetheless, the limited 
overlap between shallow and deep communities indicates 
that deep reef may still play an important role in ensuring 
local protection of certain genotypes (e.g. the small num- 
ber of A agaricites colonies that occur in the deep) when 
adjacent shallow- reefs get affected by a major disturbance. 

Conclusions 

The findings in this study highlight that depth-related 
parameters have played a major role in the diversifica- 
tion of brooding corals and their associated Symbiodi- 
nium, and corroborate the notion that brooding corals 
and their Symbiodinium are engaged in processes of 
coevolution. Further studies should investigate whether 
any adaptive differentiation of these species has occurred 
at the population level, and whether the high levels of 
specificity observed here are also prevalent in other 
brooding coral genera. Overall, we need a better under- 
standing of the actual mechanisms through which the 
depth-related environmental gradients have led to popu- 
lation/species divergence on coral reefs and how this 
evolutionary history affects the current day and future 
interconnectedness of shallow and deep reef systems. 

Additional file 



Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

PB designed and conceived of the study, collected specimens, carried out 
labwork, performed the genetic analyses, and drafted the manuscript. PRF 
designed and conceived of the study, collected specimens, carried out 
labwork, helped with the interpretation of data, and drafted the manuscript. 
JJO, KBH and NE collected specimens, carried out labwork and helped to 
draft the manuscript. JvB carried out labwork, and helped to draft the 
manuscript. MJHV carried out the transects, collected specimens and helped 
to draft the manuscript. RPMB, PMV, and OHG designed and conceived of 
the study and helped to draft the manuscript. All authors read and approved 
the final manuscript. 



Acknowledgements 

The authors thank Floris de Jongh, John Crocker, Dennis Hamro-Drotz, and 
other students that have helped over the years with the collection of coral 
specimens. We thank Judith Bakker and Joshua Boldt for their support 
during deep dives, Carlos "No joh" Winterdaal for tremendous support in the 
field, Harry Witte and Linda Tonk for their assistance in the lab, and Maria 
Bongaerts for providing input on manuscript figures. This study was funded 
by the Australian Research Council (ARC) Centre of Excellence for Coral Reef 
Studies, Pacific Blue Foundation, Schure-Beijerinck-Popping fund (KNAW), 
and supported by the CARMABI Research Station. 

Author details 

School of Biological Sciences, The University of Queensland, 4072 St Lucia, 
QLD, Australia. 2 Global Change Institute, The University of Queensland, 4072 
St Lucia, QLD, Australia. 3 ARC Centre of Excellence for Coral Reef Studies, The 
University of Queensland, 4072 St Lucia, QLD, Australia. 4 CARMABI, 
Piscaderabaai z/n, PO Box 2090, Willemstad, Curagao. department of Marine 
Biology, University of Vienna, Althanstrasse 14, 1090 Vienna, Austria. 
Netherlands Institute for Sea Research, PO Box 59, 1790 Den Burg, AB, The 
Netherlands. 7 Heron Island Research Station, The University of Queensland, 
4680 Heron Island, QLD, Australia. 8 Aquatic Microbiology, Institute for 
Biodiversity and Ecosystem Dynamics, University of Amsterdam, PO Box 
94248, 1090 Amsterdam, GE, The Netherlands. 

Received: 15 June 2012 Accepted: 14 August 2013 
Published: 23 September 2013 

References 

1. Connell JH: Diversity in tropical rain forests and coral reefs. Science 1978, 

199:1302-1310. 

2. Knowlton N, Jackson JBC: New taxonomy and niche partitioning on coral 
reefs: jack of all trades or master of some? Trends Ecol Evol 1994, 1:7-9. 

3. Souter P: Hidden genetic diversity in a key model species of coral. 
Mar Biol 2010, 157:875-885. 

4. Stefani F, Benzoni F, Yang S-Y, Pichon M, Galli P, Chen CA: Comparison of 
morphological and genetic analyses reveal cryptic divergence and 
morphological plasticity in Stylophora (Cnidaria, Scleractinia). Coral Reefs 
2011,30:1033-1049. 

5. Flot J-F, Blanchot J, Charpy L, Cruaud C, Licuanan WY, Nakano Y, Payri C, 
Tillier S: Incongruence between morphotypes and genetically delimited 
species in the coral genus Stylophora: phenotypic plasticity, 
morphological convergence, morphological stasis or interspecific 
hybridization? BMC Ecol 201 1 , 1 1 :22. 

6. Carlon DB, Budd AF: Incipient speciation across a depth gradient in a 
scleractinian coral? Evolution 2002, 56:2227-2242. 

7. Vermeij MJA, Sandin SA, Samhouri JF: Local habitat composition 
determines the relative frequency and interbreeding potential for two 
Caribbean coral morphospecies. Evol Ecol 2007, 21:27-47. 

8. Bongaerts P, Riginos C, Ridgway T, Sampayo EM, Van Oppen MJH, Englebert 
N, Vermeulen F, Hoegh-Guldberg O: Genetic divergence across habitats in 
the widespread coral Seriatopora hystrix and its associated 
Symbiodinium. PLoS One 2010, 5 :e 10871. 

9. Frade PR, Reyes-Nivia MC, Faria J, Kaandorp JA, Luttikhuizen PC, Bak RPM: Semi- 
permeable species boundaries in the coral genus Madracis: introgression in 
a brooding coral system. Mol Phylogen Evo/ 2010, 57(3):1 072-1 090. 

10. Todd PA: Morphological plasticity in scleractinian corals. Biol Rev 2008, 
83:315-337. 

11. Ow XY, Todd PA: Light-induced morphological plasticity in the 
scleractinian coral Goniastrea pectinata and its functional significance. 
Coral Reefs 20]0, 29:797-808. 

12. Freudenthal HD: Symbiodinium gen. nov. and Symbiodinium 
microadriaticum sp. nov, a zooxanthella: taxonomy, life cycle and 
morphology. J Protozool 1962, 9:45-53. 

13. Rowan R, Powers DA: A molecular genetic classification of zooxanthellae 
and the evolution of animal-algal symbiosis. Science 1991, 251:1348-1351. 

14. LaJeunesse TC: Investigating the biodiversity, ecology, and phylogeny of 
endosymbiotic dinoflagellates in the genus Symbiodinium using the ITS 
region: in search of a 'species' level marker. J Phycol 2001, 37(5):866-880. 

15. Takabayashi M, Santos SR, Cook CB: Mitochondrial DNA phylogeny of the 
symbiotic dinoflagellates {Symbiodinium, Dinophyta). J Phycol 2004, 
40(1 ):1 60-1 64. 



Additional file 1: Denaturing gradient gel electrophoresis of 
Symbiodinium ITS2 types associated with Agaricia species. 

Sequences used to characterize each symbiont profile are shown 
adjacent to bands in the gel image. Types in italics represent novel 
sequences, with the name specifying the sequence to which they are 
most related, followed by a capital N (indicating novel sequence) and an 
arbitrary number (e.g. C3b/V7). 



Bongaerts et al. BMC Evolutionary Biology 2013, 13:205 
http://www.biomedcentral.eom/1 47 1 -2 1 48/1 3/205 



16. LaJeunesse TC: "Species" radiations of symbiotic dinoflagellates in the 
Atlantic and Indo-Pacific since the Miocene-Pliocene transition. 

Mol Biol Evol 2005, 22:570-581. 

17. LaJeunesse TC, Parkinson JE, Reimer JD: A genetics-based description of 
Symbiodinium minutum sp. nov. and S. psygmophilum sp. nov. 
(dinophyceae), two dinoflagellates symbiotic with cnidaria. J Phycol 2012, 
48:1380-1391. 

18. Ulstrup KE, Van Oppen MJ: Geographic and habitat partitioning of 
genetically distinct zooxanthellae {Symbiodinium) in Acropora corals on 
the Great Barrier Reef. Mol Ecol 2003, 12:3477-3484. 

19. Iglesias-Prieto R, Beltran VH, LaJeunesse TC, Reyes-Bonilla H, Thome PE: 
Different algal symbionts explain the vertical distribution of dominant 
reef corals in the eastern Pacific. Proc R Soc B 2004, 271:1757-1763. 

20. Sampayo EM, Franceschinis L, Hoegh-Guldberg O, Dove S: Niche 
partitioning of closely related symbiotic dinoflagellates. Mol Ecol 2007, 
16:3721-3733. 

21. Frade PR, De Jongh F, Vermeulen F, Van Bleijswijk J, Bak RPM: Variation in 
symbiont distribution between closely related coral species over large 
depth ranges. Mol Ecol 2008, 17:691-703. 

22. LaJeunesse TC, Pettay DT, Sampayo EM, Phongsuwan N, Brown B, 
Obura DO, Hoegh-Guldberg O, Fitt WK: Long-standing environmental 
conditions, geographic isolation and host-symbiont specificity 
influence the relative ecological dominance and genetic 
diversification of coral endosymbionts in the genus Symbiodinium. 

J Biogeogr 2010, 37:785-800. 

23. LaJeunesse TC, Thornhill DJ, Cox EF, Stanton FG, Fitt WK, Schmidt GW: High 
diversity and host specificity observed among symbiotic dinoflagellates 
in reef coral communities from Hawaii. Coral Reefs 2004, 23:596-603. 

24. Loh WKW, Loi T, Carter D, Hoegh-Guldberg O: Genetic variability of the 
symbiotic dinoflagellates from the wide ranging coral species 
Seriotopora hystrix and Acropora longicyothus in the Indo-West Pacific. 
Mar Ecol Prog Ser 2001, 222:97-107. 

25. Stat M, Hoegh-Guldberg O, Fitt WK, Carter D: Host symbiont acquisition 
strategy drives Symbiodinium diversity in the southern Great Barrier 
Reef. Coral Reefs 2008, 27:763-772. 

26. Bongaerts P, Riginos C, Hay K, Van Oppen MJH, Hoegh-Guldberg O, Dove S: 
Adaptive divergence in a scleractinian coral: physiological adaptation of 
Seriatopora hystrix to shallow and deep reef habitats. BMC Evol Biol 201 1, 
11:303. 

27. Bongaerts P, Sampayo EM, Bridge TLC, Ridgway T, Vermeulen F, Englebert 
N, Webster JM, Hoegh-Guldberg O: Symbiodinium diversity of mesophotic 
coral communities on the GBR: a first assessment. Mar Ecol Prog Ser 201 1, 
439:117-126. 

28. LaJeunesse TC, Smith RT, Finney J, Oxenford H: Outbreak and persistence 
of opportunistic symbiotic dinoflagellates during the 2005 Caribbean 
mass coral 'bleaching' event. Proc R Soc B 2009, 276:4139-4148. 

29. Van Oppen MJH: Mode of zooxanthella transmission does not affect 
zooxanthella diversity in acroporid corals. Mar Biol 2004, 144:1-7. 

30. Hoegh-Guldberg O, Mumby PJ, Hooten AJ, Steneck RS, Greenfield P, Gomez E, 
Harvell CD, Sale PF, Edwards AJ, Caldeira K, Knowlton N, CEakin M, Iglesias-Prieto R, 
Muthiga N, Bradbury RH, Dubi A, Hatziolos ME: Coral reefs under rapid climate 
change and ocean acidification. Science 2007, 318:1 737-1 742. 

31. Bongaerts P, Ridgway T, Sampayo EM, Hoegh-Guldberg O: Assessing the 
'deep reef refugia' hypothesis: focus on Caribbean reefs. Coral Reefs 2010, 
29(2):309-327. 

32. Van Oppen MJH, Bongaerts P, Underwood J, Peplow L, Cooper T: The role 
of deep reefs in shallow reef recovery: an assessment of vertical 
connectivity in a brooding coral from west and east Australia. 

Mol Ecol 2011,20:1647-1660. 

33. Rowan R, Knowlton N: Intraspecific diversity and ecological zonation in 
coral-algal symbiosis. Proc Natl Acad Sci USA 1995, 92:2850-2853. 

34. Frade PR, Bongaerts P, Winkelhagen AJS, Tonk L, Bak RPM: In situ 
photobiology of corals over large depth ranges: a multivariate analysis 
on the roles of environment, host, and algal symbiont. Limnol Oceanogr 
2008,53:2711-2723. 

35. Lesser MP, Slattery M, Stat M, Ojimi M, Gates RD, Grottoli A: 
Photoacclimatization by the coral Montastraea cavernosa in the 
mesophotic zone: light, food, and genetics. Ecology 201 0, 91 :990-1 003. 

36. Richmond RH, Hunter CL: Reproduction and recruitment of corals: 
Comparisons among the Caribbean, the Tropical Pacific, and the Red 
Sea. Mar Ecol Prog Ser 1990, 60:185-203. 



Page 13 of 14 



37. Goreau TF, Goreau Nl: The ecology of Jamaican coral reefs. II. 
Geomorphology, zonation, and sedimentary phases. Bull Mar Sci 1973, 
23:399-464. 

38. Hughes TP, Jackson JBC: Population dynamics and life histories of 
foliaceous corals. Ecol Monogr 1985, 55:141-166. 

39. Bak RPM, Nieuwland G, Meesters EH: Coral reef crisis in deep and shallow 
reefs: 30 years of constancy and change in reefs of Curacao and Bonaire. 
Coral Reefs 2005, 24:475-479. 

40. Kahng SE, Garcia R, Spalding HL, Brokovich E, Wagner D, Weil E, Hinderstein 
L, Toonen RJ: Community ecology of mesophotic coral reef ecosystems. 
Coral Reefs 2010, 29:255-275. 

41. Titlyanov EA: Structure and morphological differences of colonies of reef- 
building branched corals from habitats with different light conditions. 
Mar Biol 1987, 1:32-36. 

42. Van Moorsel GWNM: Reproductive strategies in two closely related stony 
corals {Agaricia, Scleractinia). Mar Ecol Prog Ser 1983, 13:273-283. 

43. Bak RPM: The growth of coral colonies and the importance of crustose 
coralline algae and burrowing sponges in relation with carbonate 
accumulation. Neth J Sea Res 1976, 10:285-337. 

44. Bak RPM, Engel MS: Distribution, abundance and survival of juvenile 
hermatypic corals (Scleractinia) and the importance of life history 
strategies in the parent coral community. Mar Biol 1979, 54:341-352. 

45. Vermeij MJ, Bakker J, Hal N, Bak RP: Juvenile coral abundance has 
decreased by more than 50% in only three decades on a small 
Caribbean island. Diversity 201 1, 3(3):296-307. 

46. LaJeunesse T: Diversity and community structure of symbiotic 
dinoflagellates from Caribbean coral reefs. Mar Biol 2002, 141:387-400. 

47. Warner ME, LaJeunesse TC, Robison JD, Thur RM: The ecological 
distribution and comparative photobiology of symbiotic dinoflagellates 
from reef corals in Belize: potential implications for coral bleaching. 
Limnol Oceanogr 2006, 51:1887-1897. 

48. Banaszak AT, Santos MGB, LaJeunesse TC, Lesser MP: The distribution of 
mycosporine-like amino acids (MAAs) and the phylogenetic identity of 
symbiotic dinoflagellates in cnidarian hosts from the Mexican Caribbean. 
J Exp Mar Bio Ecol 2006, 337(2):1 31 -1 46. 

49. Santos RS, LaJeunesse TC: Searchable database of Symbiodinium diversity- 
geographic and ecological diversity (SD2-GED); 2006. http://www.auburn.edu/ 
~santosr/sd2_ged.htm. 

50. Finney JC, Pettay DT, Sampayo EM, Warner ME, Oxenford HA, LaJeunesse 
TC: The relative significance of host-habitat, depth, and geography on 
the ecology, endemism, and speciation of coral endosymbionts in the 
genus symbiodinium. Microb Ecol 2010, 60:250-263. 

51. Frade PR, Englebert N, Faria J, Visser PM, Bak RPM: Distribution and 
photobiology of Symbiodinium types in different light environments for 
three colour morphs of the coral Madracis pharensis: is there more to it 
than total irradiance? Coral Reefs 2008, 27:913-925. 

52. Bak RPM: Coral reefs and their zonation in Netherlands Antilles. Stud Geol 
1977, 4:3-16. 

53. Wells JW: New and old scleractinian corals from Jamaica. Bull Mar Sci 
1973, 23:16-54. 

54. Veron JEN, Stafford-Smith M: Corals of the World. Townsville: Australian 
Institute of Marine Science; 2000. 

55. Humann P, DeLoach N: Reef coral identification. New World, Jacksonville, 
Florida: Florida Caribbean Bahamas including marine plants; 2002. 

56. Wilson K, Li Y, Whan V, Lehnert S, Byrne K, Moore S, Pongsomboon S, 
Tassanakajon A, Rosenberg G, Bailment E: Genetic mapping of the black 
tiger shrimp Penaeus monodon with amplified fragment length 
polymorphism. Aquaculture 2002, 204:297-309. 

57. Mieog JC, Van Oppen MJH, Cantin NE, Stam WT, Olsen JL: Real-time PCR 
reveals a high incidence of Symbiodinium clade D at low levels in four 
scleractinian corals across the Great Barrier Reef: implications for 
symbiont shuffling. Coral Reefs 2007, 26:449-457. 

58. Silverstein RN, Correa AMS, Baker AC: Specificity is rarely absolute in 
coral-algal symbiosis: implications for coral response to climate change. 
Proc Roy Soc B 2012. doi:1 0.1 098/rspb.201 2.0055. 

59. Tamura K, Dudley J, Nei M, Kumar S: MEGA4: molecular evolutionary 
genetics analysis (MEGA) software version 4.0. Mol Biol Evol 2007, 
24:1596-1599. 

60. Medina M, Collins AG, Takaoka TL, Kuehl JV, Boore JL: Naked corals: 
skeleton loss in Scleractinia. Proc Natl Acad Sci USA 2006, 
103(24):9096-9100. 



Bongaerts et al. BMC Evolutionary Biology 201 3, 1 3:205 Page 1 4 of 1 4 

http://www.biomedcentral.eom/1 47 1 -2 1 48/1 3/205 



61. 



62. 



63. 



64. 



65. 



66. 



67. 



69. 



70. 



71. 



72. 



73. 



74. 



75. 



76. 



77. 



79. 



Shearer TL, Coffroth MA: DNA Barcoding: Barcoding corals: limited by 
interspecific divergence, not intraspecific variation. Mol Ecol Resour 2008, 
8(2):247-255. 

Fukami H, Chen CA, Budd AF, Collins A, Wallace C, Chuang YY, Chen C, Dai 
CF, Iwao K, Sheppard C, Knowlton N: Mitochondrial and nuclear genes 
suggest that stony corals are monophyletic but most families of stony 
corals are not (Order Scleractinia, Class Anthozoa, Phylum Cnidaria). 

PLoS One 2008, 3:e3222. 

Barbeitos MS, Romano SL, Lasker HR: Repeated loss of coloniality and 
symbiosis in scleractinian corals. Proc Natl Acad Sci USA 2010, 
107(26):1 1877-1 1882. 

Posada D: jModelTest: phylogenetic model averaging. Mol Biol Evol 2008, 
25(7):1 253-1 256. 

Peakall R, Smouse PE: GenAlEx V5: Genetic Analysis in Excel. Population 
Genetic Software for Teaching and Research: Australian National University; 
2001. 

Thornhill DJ, LaJeunesse TC, Santos SR: Measuring rDNA diversity in 
eukaryotic microbial systems: how intragenomic variation, pseudogenes, 
and PCR artifacts confound biodiversity estimates. Mol Ecol 2007, 
16:5326-5340. 

LaJeunesse TC, Smith R, Walther M, Pinzon J, Pettay DT, McGinley M, 
Aschaffenburg M, Medina-Rosas P, Cupul-Magaha AL, Lopez Perez A, 
Reyes-Bonilla H, Warner ME: Host- symbiont recombination versus 
natural selection in the response of coral - dinoflagellate symbioses to 
environmental disturbance. Proc Roy Soc B 2010, 277:2925-2934. 
Hoeksema BW: Evolutionary trends in onshore-offshore distribution 
patterns of mushroom coral species (Scleractinia: Fungiidae). Contrib Zool 
2012,81:199-221. 

Goldberg WM: The ecology of the coral octocoral communities off the 
southeast Florida coast: geomorphology, species composition and 
zonation. Bull Mar Sci 1 973, 23:465-488. 

Van den Hoek C, Breeman AM, Bak RPM, Van Buurt G: The distribution of 

algae, corals, and gorgonians in relation to depth, light attenuation, 

water movement, and grazing pressure in the fringing coral reef of 

Curacao, Netherlands Antilles. Aquat Bot 1978, 5:1-46. 

Bak RPM, Luckhurst BE: Constancy and change in coral reef habitats along 

depth gradients at Curacao. Oecologia 1980, 47:145-155. 

Kuhlmann D: Composition and ecology of deep-water coral associations. 

Helgol Mar Res 1983, 36:183-204. 

Chan Y, Pochon X, Fisher MA, Wagner D, Concepcion GT, Kahng SE, Toonen 
RJ, Gates RD: Generalist dinoflagellate endosymbionts and host genotype 
diversity detected from mesophotic (67-100 m depths) coral Leptoseris. 

BMC Ecol 2009, 9:21. 

Mass T, Kline Dl, Roopin M, Veal CJ, Cohen S, Muz D, Levy O: The spectral 
quality of light is a key driver of photosynthesis and photoadaptation in 
Stylophora pistillata colonies from different depths in the Red Sea. 

J Exp Biol 2010, 213:4084-4091. 

Winters G, Beer S, Zvi BB, Brickner I, Loya Y: Spatial and temporal 
photoacclimation of Stylophora pistillata: zooxanthella size, 
pigmentation, location and clade. Mar Ecol Prog Ser 2009, 384:107-1 19. 
Vermeij MJA, Bak RPM: Species-specific population structure of closely 
related coral morphospecies along a depth gradient (5-60 m) over a 
Caribbean reef slope. Bull Mar Sci 2003, 73:725-744. 
Van Oppen MJ, Palstra FP, Piquet AM, Miller DJ: Patterns of coral- 
dinoflagellate associations in Acropora: significance of local availability 
and physiology of Symbiodinium strains and host-symbiont selectivity. 
Proc R Soc B 2001, 268:1759-1767. 

Van Oppen MJ, Willis BL, Van Rheede T, Miller DJ: Spawning times, 
reproductive compatibilities and genetic structuring in the Acropora 
aspera group: evidence for natural hybridization and semi-permeable 
species boundaries in corals. Mol Ecol 2002, 1 1:1363-1376. 
Fukami H, Budd AF, Levitan DR, Jara J, Kersanach R, Knowlton N: 
Geographic differences in species boundaries among members of the 
Montastraea annularis complex based on molecular and morphological 
markers. Evolution 2004, 58:324-337. 

Flot JF, Licuanan WY, Nakano Y, Payri C, Cruaud C, Tillier S: Mitochondrial 
sequences of Seriatopora corals show little agreement with morphology 
and reveal the duplication of a tRNA gene near the control region. 

Coral Reefs 2008, 27:789-794. 



81. Forsman ZH, Barshis DJ, Hunter CL, Toonen RJ: Shape-shifting corals: 
molecular markers show morphology is evolutionarily plastic in Pontes. 
BMC Evol Biol 2009, 9:45. 

82. Miller KJ, Ayre DJ: The role of sexual and asexual reproduction in 
structuring high latitude populations of the reef coral Pocillopora 
damicornis. Heredity 2004, 92:557-568. 

83. Diekmann O, Olsen J, Stam W, Bak R: Genetic variation within 
Symbiodinium clade B from the coral genus Madracis in the Caribbean 
(Netherlands Antilles). Coral Reefs 2003, 22:29-33. 

84. Diekmann OE, Bak RPM, Stam WT, Olsen JL: Molecular genetic evidence 
for probable reticulate speciation in the coral genus Madracis from a 
Caribbean fringing reef slope. Mar Biol 2001, 139:221-233. 

85. Vermeij MJA, Bak RPM: How are coral populations structured by light? 
Marine light regimes and the distribution of Madracis. Mar Ecol Prog Ser 
2002, 233:105-116. 

86. Vermeij MJA, Sampayo E, Broker K, Bak RPM: Variation in planulae release 
of closely related coral species. Mar Ecol Prog Ser 2003, 247:75-84. 

87. Vermeij MJA, Sampayo E, Broker K, Bak RPM: The reproductive biology of 
closely related coral species: gametogenesis in Madracis from the 
southern Caribbean. Coral Reefs 2004, 23:206-214. 

88. Ayre DJ, Hughes TP: Genotypic diversity and gene flow in brooding and 
spawning corals along the Great Barrier Reef, Australia. Evolution 2000, 
54:1590-1605. 

89. Baums IB, Miller MW, Hellberg ME: Regionally isolated populations of 
an imperiled Caribbean coral, Acropora palmata. Mol Ecol 2005, 
14:1377-1390. 

90. Underwood JN, Smith LD, Van Oppen MJH, Gilmour JP: Multiple scales of 
genetic connectivity in a brooding coral on isolated reefs following 
catastrophic bleaching. Mol Ecol 2007, 16:771-784. 



doi:1 0.1 1 86/1 471 -21 48-1 3-205 

Cite this article as: Bongaerts et al: Sharing the slope: depth 
partitioning of agariciid corals and associated Symbiodinium across 
shallow and mesophotic habitats (2-60 m) on a Caribbean reef. BMC 

Evolutionary Biology 201 3 1 3:205. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



